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■ Abstract 
^ \ 

, An additional Z§ symmetry hidden in the fermion and Higgs sec- 

^ ' tors of the Standard Model has been found recently pQ. A lattice reg- 

ularization of the Standard Model was constructed that possesses this 

\ symmetry. In j2] we have reported our results on the numerical simula- 

\^ | tion of the electroweak sector of the model. In this paper we report our 

■ results on the numerical simulation of the full (SU (3) (g) SU(2) <g> U (1) 
, ) model. The phase diagram of the model has been investigated using 

^ static quark and lepton potentials. Various types of monopoles have 

■ been constructed. Their densities appear to be sensitive to the phase 
transition lines. Differences between the realizations of the Standard 
Model which do or do not possess the mentioned Zg symmetry, are 

'"" J" 1 1 discussed. 

a: 

<D . 
43 ■ 



1 Introduction 

Until recently it was thought that all the symmetries of the Standard Model 
(SM), which must be used when dealing with its discretization, are known. 
However, in 1] it was shown that there exists an additional Z 6 = Z 2 ® Z 3 
symmetry in the fermion and Higgs sectors of the SM. It is connected to 
the centers Z 3 and Z 2 of the SU(3) and SU{2) subgroups 1 . The gauge 
sector of the SM (in its discretized form) was redefined in such a way that 
it has the same naive continuum limit as the original one, while keeping the 
mentioned symmetry. The resulting model differs from the conventional SM 
via its symmetry properties. Therefore we expect, that nonperturbatively 
these two models may represent different physics. 

1 The emergence of Z§ symmetry in the SM and its supersymmetric extension was 
independently considered in a different context in 
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Investigation of the electroweak sector of the SM with the additional 
Zq symmetry shows, that there are indeed certain differences between this 
discretization and the conventional one [2|. Namely, it has been found that 
the phase transition lines corresponding to the U(l) and SU(2) degrees of 
freedom join in a triple point, forming a common line. In contrast to this, 
in the conventional model the phase transition line corresponding to 577(2) 
degrees of freedom has an endpoint and the transition becomes continuous in 
a certain region of coupling constants [3] . In this paper we report our results 
on the full SM (including SU(3) degrees of freedom) and claim that the same 
phenomenon takes place here. Now the 577(3), SU{2) and £7(1) degrees of 
freedom are connected via their centers. This, in our opinion, is the reason 
why the phase transition lines corresponding to the phase transitions in pure 
U(l) and SU(2) models again join together forming a common line. It turns 
out that SU(3) fields experience this common phase transition as well. 

This paper is organized as follows. In the next section we summarize 
the formulation of the SM in terms of link variables and demonstrate the 
emergence of an additional Z e symmetry in its fermion and Higgs sectors. In 
Sect. El we detail the model with explicit Z 6 symmetry on the lattice, while 
in Sect.|3]we recall the definition of the maximal center projection. The next 
section contains the definitions of the quantities we measure on the lattice; 
it is followed by Sect. |H1 where we show our numerical results. We end with 
a summary. 

2 Zq symmetry in the Standard Model 

In this section we remind the reader of what we call the additional Z 6 sym- 
metry. The SM contains the following variables: 
1. The gauge field U = (T, U, 8), where 



3. Anticommuting spinor variables, representing leptons and quarks: 



r G S77(3), U E SU{2), e w e 17(1), 



(1) 



realized as link variables on the lattice. 
2. A scalar doublet 

$ a , a = 1,2. 



(2) 
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left-handed leptons 




right-handed leptons 


TUe¥ 


left-handed quarks 


1 e 3 ° 


right-handed d, s, and, b - quarks 


T-i —8 

1 es 


right-handed u, c, and, t - quarks 


Ue li> 


the Higgs scalar field 



The action has the form 

S = S g + S H + S f , (4) 

where we denote the fermion part of the action by Sf, the pure gauge part 
is denoted by S g , and the scalar part of the action by Sa- 
in any lattice realization of Sh and Sf both these terms depend upon 
link variables U considered in the representations corresponding to quarks, 
leptons, and the Higgs scalar field, respectively. Therefore U appears in 
the combinations shown in the table. Our observation is that all the listed 
combinations are invariant under the following transformations: 

U -> Ue- mN , 
9 -> 6 + nN, 

r _> re^' l ^ N , (5) 

where N is an arbitrary integer link variable. It represents a three-dimensional 
hypersurface on the dual lattice. Both Sh and Sf (in any realization) are in- 
variant under the simultaneous transformations This symmetry reveals 
the correspondence between the centers of the SU(2) and SU(3) subgroups 
of the gauge group. 

After integrating out fermion and scalar degrees of freedom any physical 
variable should depend upon gauge invariant quantities only. Those are the 
Wilson loops: UJ SU (3){C) = Trn UnkgC r link , uj SU (2)(C) = Trrii inkeC £7i ink , and 
u>u(i)(C) = riii nkgC exp(|6 l ii nk ). Here C is an arbitrary closed contour on the 
lattice (with self - intersections allowed). These Wilson loops are trivially 
invariant under the transformation (0) with the field N representing a closed 
three-dimensional hypersurface on the dual lattice. Therefore the nontrivial 
part of the symmetry (J3J) corresponds to a closed two-dimensional surface 
on the dual lattice that is the boundary of the hypersurface represented by 



3 



N. Then in terms of the gauge invariant quantities u the transformation © 
acquires the form: 

Uu(i)(C) -> exp(-i±vrL(C, Y,)) uj u{1) (C) 
usu{2){C) -> exp(i7rL(C,S)) ^5(7(2) (C) 
^5c/(3)(C) -> exp(z§7rL(C,E))u;s[/ ( 3)(C) (6) 

Here £ is an arbitrary closed surface (on the dual lattice) and L(C, E) is the 
integer linking number of this surface and the closed contour C. From © it 
follows, that the symmetry is of Z e type. 



3 The model under investigation 

It is obvious that the pure gauge-field part of the action in its conventional 
continuum formulation (or, say, in lattice Wilson formulation) is not invariant 
under (JBJ). However, the lattice realization of the pure gauge field term of the 
action can be constructed in such a way that it also preserves the mentioned 
symmetry. For the reasons listed in pQ we consider it in the following form: 

Sg = P E {2(l-fTr[/ p cos0 p ) + 

plaquettes 

+ (1 -cos2# p ) 

+ 6[1 - ±ReTrr p Tr[/ p exp(^ p /3)] 

+ 3[1 - ±ReTrr p exp(-2i0 p /3)] 

+ 3[1 - |ReTrr p exp(4^ p /3)]}, (7) 

where the sum runs over the elementary plaquettes of the lattice. Each 
term of the action Eq. J7J) corresponds to a parallel transporter along the 
boundary dp of plaquette p. The corresponding plaquette variables con- 
structed of lattice gauge fields are U p = u)su{2)(dp) ,T P = ujsu(3)(9p) , and 
6 P = Argu u{1) (dp). 

The potential for the scalar field is considered in its simplest form [2] in 
the London limit, i.e., in the limit of infinite bare Higgs mass. After fixing 
the unitary gauge we obtain: 

S H = 7 E[l-&(^)]. (8) 

xy 
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The following variables are (naively) considered as creating a photon, Z 
boson, and W boson respectively: 

A xy = A% = [AigU^ y +6 xy }mod2n, 
Z xy = Z» = [AigU£-e xu ]mod2n, 
W xy = W£ = Ul 2 y e w *y. (9) 

Here, \i represents the direction (xy). After fixing the unitary gauge the 
electromagnetic U(l) symmetry remains: 

U xy > 9 x U X y9yi 

xy -> 9 xy + a y /2-a x /2, (10) 
where g x = diag(e tax ^ 2 , e~ l0tx l 2 ). The fields A, Z, and W transform as follows: 



A 


~^ A xy ~\~ oc y (^xi 




7 

^-*xy 


t-'xyi 




w xy - 


-> W xy e~ ia \ 


(11) 



We consider our model in quenched approximation, i.e., we neglect the 
effect of virtual fermion loops. Therefore the particular form of Sf is not of 
interest for us at this stage. 

In order to extract physical information from the 577(3) fields in a par- 
ticulary simple way we use the so-called indirect Maximal Center Projection 
(see, for example, 0E])- 



4 The Maximal Center Projection. 

The Maximal Center Projection makes the link matrix T as close as possible 
to the elements of the center Z 3 of SU(3): Z 3 = {diag(e^ 2wi ^ N , e (2 ™M N , e^ 3 ^}, 
where N 6 {1, 0, —1}. The procedure works as follows. 
First, make the functional 

Qi = E(l r nl + l r 22l + |r 3 3|) (12) 

links 

maximal with respect to the gauge transformations T xy — ► g].r xy g y , thus fix- 
ing the Maximal Abelian gauge. As a consequence every link matrix becomes 
almost diagonal. 
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Secondly, to make this matrix as close as possible to the center of 577(3), 
make the phases of the diagonal elements of this matrix maximally close to 
each other. This is done by minimizing the functional 

Q 2 = ^{[l-cos(Arg(r 11 )-Arg(r 22 ))] + [l-cos(Arg(r 11 )-Arg(r 3 3))] 

links 

+ [1 - cos(Arg(r 22 ) - Arg(r 33 ))]}. (13) 

with respect to the gauge transformations. This gauge condition is invariant 
under the central subgroup Z 3 of 577(3). 

In our model 577(3) fields are connected with the £7(1) and SU(2) fields 
via the center of the gauge group. Therefore, instead of the center vortices 
and center monopoles we define various kinds of monopole - like fields. The 
definitions of these fields includes the following integer-valued link variable 
N (defined after fixing the Maximal Center gauge): 

N xy = if (Arg(r 11 )+Arg(r 22 )+Arg(r 33 ))/3e]-7r/3,7r/3], 
N xy = 1 if (Arg(r n ) + Arg(r 22 ) + Arg(r 33 ))/3 G ]vr/3, tt], 
A^ = -l if (Arg(r 11 )+Arg(r 22 )+Arg(r 33 ))/3e]-7r,-7r/3]. (14) 

In other words, N = if T is close to 1, N = 1 if T is close to e 2nl ^ 3 and 
N = -1 if T is close to e" 2 ™/ 3 . 

Next, we define the following link fields 



r 1 - 

xy 


g ™xy 


+ AxgUH + -9 xy mod27r, 




r 2 - 

^ xy 




r^ u xy 


mod27r, 




C xy = 


^A 
. 3 xy 


+ -®xy 


mod27r. 


(15) 



These fields correspond to the last three terms of Eq. ((7|). Their construction 
comes from the representation of Y as a product of exp((27ri/3)A) and V, 
where V is the SU(3)/Z 3 variable (Arg(Vn) + Arg(V" 22 ) + Arg(^ 33 ))/3 G 
] — 7r/3,7r/3]. Thus T = exp((2iii/3)N)V. We expect, that ((7j) suppresses 
Vpi aq and Cp laq , i = 1,2,3 while the fields N, |, and U\\ (being considered 
independently of each other) are expected to be disordered. This assumption 
is justified by the numerical simulations. 
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5 Quantities to be measured 

We investigated five types of monopoles. The monopoles, which carry infor- 
mation about colored fields are extracted from C % : 

j ci = — *d([dC i ]mod2ir). (16) 
2n 

(Here we used the notations of differential forms on the lattice. For a defini- 
tion of those notations see, for example, 

Pure U(l) monopoles, corresponding to the second term in (JJJ), are ex- 
tracted from 29: 

he = ^-*d{[d26]mod2n). (17) 
The electromagnetic monopoles, corresponding to the first term in (jZJ), are: 

J A = ^*d([dA]mod2n). (18) 
The density of the monopoles is defined as follows: 

/ Slinks Jlink \ /-i n\ 

P= \ 4L* /' ^ 

where L is the lattice size. To understand the dynamics of external charged 
particles, we consider the Wilson loops defined in the fermion representations 
listed above (in the table): 

WSptCO = (ReTTU (xy)el U xy e-^y), 
W£ pt (0 = (ReU {xy)el e- 2 ^y), 
W q L uarks (/) = (ReU {xyW T xy U xy e^), 

W d R own quarks (0 = (^U {xy)el ? xy e"T »»» ) , 

W u R pquarks (0 = (Ren (xy)el T xy e^). (20) 

Here / denotes a closed contour on the lattice. We consider the following 
quantity constructed from the rectangular Wilson loop of size a x t: 

V(a) = hm , W (°, X *> -. (21) 
v ; t^oo W(a x (t + 1)) V ; 

A linear behavior of V(a) would indicate the existence of a charge - anti 
charge string with nonzero tension. 
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6 Numerical results 



In our calculations we investigated lattices L 4 for L = 6, L = 12, and L = 16 
with symmetric boundary conditions. 

We summarize our qualitative results in the phase diagram represented 
in Fig. ^ The model contains three phases. The first one (I) is a phase, in 
which the dynamics of external leptons is confinement - like, i.e. is similar 
to that of external charges in QCD with dynamical fermions. In the second 
phase (II) the behavior of left-handed leptons is confinement-like, while for 
right-handed ones it is not. The last one (III) is the Higgs phase, in which no 
confining forces between leptons are observed at all. In all three phases there 
is the confinement of all external quark fields (left quarks, right up quarks, 
right down quarks). 

This is illustrated by Figs. El in which we show V(a) extracted from 
the Wilson loops Eq. (|2Uj) at two typical points that belong to phases II 
(7 = 0.5) and III (7 = 1.5) of the model (the behavior of all potentials in the 
phase I is confinement - like). We represent here the potential for only one 
colored Wilson loop, i.e. for Wq Uarks , because the string tension extracted 
from the other two potentials coincides with the string tension extracted 
from the potential represented in the figure within the errors. This is, of 
course, exactly what we have expected: string tensions for different types 
of quarks are equal to each other. Thus, the potential, extracted from the 
colored fields, possesses linear behavior in all phases, indicating appearance 
of confinement of quarks. 

By making a linear fit to the lepton potentials at values a > 5 we found 
that only in the case of left-handed leptons the value of the string tension 
is much larger than its statistical uncertainty in phase II. For left-handed 
leptons in the Higgs phase and right-handed leptons in both phases, the un- 
certainty in the values of the string tension turns out to be larger than about 
24% of its value. In these cases we do not consider the string tension to be sig- 
nificantly different from zero. However, as for QCD with dynamical fermions 
or the SU(2) fundamental Higgs model [TUIIEI], these results do not mean 
that confinement of leptons occurs. The charge - anti charge string must 
be torn by virtual charged scalar particles, which are present in the vacuum 
due to the Higgs field. Thus V(a) may be linear only at sufficiently small 
distances, while starting from some distance it must not increase, indicating 
the breaking of the string. Unfortunately the accuracy of our measurements 
does not allow us to observe this phenomenon in detail. 
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The connection between the properties of monopoles and the phase struc- 
ture of the model is illustrated by Figs. |HJ which shows the monopole density 
versus j3 at fixed 7 = 0.5. Again we represent here only one type of the three 
monopoles, which have colored origin. Namely, we consider jqi. (Behavior 
of the others is similar.) One can see, that the density of the 29 - monopoles 
as well as C 1 - monopoles falls sharply in phase II, while the electromagnetic 
monopole density does not. 

We note here, that according to our measurements the electromagnetic 
monopole density falls to zero while shifting from phases I to III. The colored 
monopoles and 29 - monopole densities fall sharply in the phase III as well. 
Thus monopoles composed of colored fields feel the phase transition, which 
are due, according to our intuition, to the £7(1) variables. This happens again 
because the Z$ symmetry binds £7(1) variables with the center of the 577(3) 
subgroup of the gauge group. 

As in [2] we mention here that the SU(2) fundamental Higgs model, has 
a similar phase structure as our model, except for the absence of the phase 
transition line between phases I and II. In the latter model it was shown 
that different phases are actually not different. This means that the phase 
transition line ends at some point and the transition between two states of 
the model becomes continuous. Thus one may expect that in our model the 
phase transition line between phases I and III ends at some point. However, 
we do not observe this for the considered values of couplings. 

In our model both phase transition lines join in a triple point, forming a 
common line. This is, evidently, the consequence of the mentioned additional 
symmetry that relates SU(2) , U(l), and SU(3) excitations. The same pic- 
ture, of course, does not emerge in the conventional SU(3) ® 577(2) ® U(l) 
gauge - Higgs model: its 577(2) ® U(l) part was investigated, for example, 
in As for the SU(3) gauge theory, it has no phase transition at finite (3 
and zero temperature at all. 

7 Conclusions 

We summarize our results as follows: 

1. We performed a numerical investigation of the quenched lattice model 
that respects the additional symmetry. 
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2. The lattice model contains three phases. In the first phase the potential 
between static leptons is confinement-like. In the second phase the 
confining forces are observed, at sufficiently small distances, between 
the left-handed external leptons. The last one is the Higgs phase, where 
there are no confinement - like forces between static leptons at all. 

3. Investigation of the monopoles constructed of colored fields shows that 
colored fields feel the phase transition lines. 

4. In all phases of the model we observe confinement of quarks. The string 
tensions for different kinds of quarks are equal to each other. 

5. The main consequence of the emergence of the additional symmetry is 
that the phase transition lines corresponding to the SU(2) and £7(1) 
degrees of freedom join in a triple point forming a common line. This 
reflects the fact that the SU(2) and U(l) excitations are related due 
to the mentioned symmetry. The same situation does not occur in the 
conventional £77(2) <S> U(l) gauge - Higgs model jlj. 

So, we have found a qualitative difference between the conventional dis- 
cretization and the discretization that respects the invariance under the 
transformations given in Eq. (JBJ). 

In order to illustrate other possible differences let us consider the problem 
of constructing the operator which creates a glueball in the Z§ invariant 
version of the lattice SM. Here we cannot use the conventional expression 

C = (l--Retrr plj J (22) 

as it is not invariant under our Z§ symmetry. Instead we may use Z§ - 
invariant expressions like 

0=(l- iRetr{r plaq e-f d v^}\ _ 1(1 _ cos2fl plaq ) (23) 

In the naive continuum limit the above expressions (|2*2*|) and (f2*3*j) coincide. In 
a similar way the naive continuum limit of the action (J7J) coincides with that 
of the conventional lattice SM action for the appropriate choice of coupling 
constants. 

However, this coincidence does not mean necessarily, that either the mod- 
els themselves or the correlators of operators (|22j) and (J23j) lead to the same 
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results. Let us recall here two precedents, i.e., two similar situations, where 
the coincidence of the naive continuum limits does not lead to the same 
physics. 

The first example is the massless lattice fermion. One may compare 
Wilson fermions with the simplest direct discretization of the Dirac fermion 
action^!]. These two actions differ from each other by a term which naively 
vanishes in the continuum limit. However, the corresponding models are not 
identical from the physical point of view. Namely, the second one contains 15 
additional fermion species while in the Wilson formulation all of them acquire 
infinite mass and disappear in the continuum limit. This phenomenon of 
fermion doubling is widely discussed in the literature. It is worth mentioning 
that another difference between these two formulations is the absence of exact 
chiral symmetry in the Wilson formulation and its appearance in the naive 
discretization. 

The second example is the pure nonabelian gauge theory. If we would 
discretize its form written in terms of gauge potentials losing the exact gauge 
invariance, the resulting lattice model would have the same naive continuum 
limit as the conventional lattice gluodynamics, which is written in terms of 
link matrices. However, in such a definition of lattice gauge theory confine- 
ment is lost [12]. 

In the two examples of lattice models considered above, which have the 
same naive continuum limit but different symmetry properties, finally lead 
to different physics. Exactly the same situation may be present in our case, 
where the naive continuum limit of the two lattice realizations of the SM is 
the same, while only one formulation is Z§ invariant. 

Another argument in favor of the point of view that these two models 
are indeed different, comes from the direct consideration of how continuum 
physics emerges in the lattice SM. Namely, there are indications [HI E] that 
several kinds of singular field configurations may survive in the continuum 
limit of nonabelian lattice gauge models. If so, the conventional action of the 
lattice SM and the action (J7J) may appear to be different when approaching 
the continuum for singular field configurations of various kinds. 
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Figure 1: The phase diagram of the model in the (f3, 7)-plane. 
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a a left-handed leptons y = 0.5 
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Figure 2: Vl(o) calculated at /3 = 0.7. Here the potentials are extracted 
from Wq Uarks (left - handed quarks), W^ pt (left - handed leptons), and W^ pt 
(right - handed leptons). 
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